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The time-dependent Hartree-Fock approach may be derived from a variational principle and a 
Slater Determinant wavefunction Ansatz. It gives a good description of nuclear processes in which 
one-body collisions dominate and has been applied with success to giant resonances and collisions 
around the barrier. It is inherently unable to give a good description of two-body observables. A 
variational principle, due to Balian and Veneroni has been proposed which can be geared to good 
reproduction of two-body observables. Keeping the Slater Determinant Ansatz, and restricting the 
two-body observables to be the squares of one-body observables, the procedure can be implemented 
as a modification of the time-dependent Hartree-Fock procedure. Applications, using the Skyrme 
effective interaction, are presented for the mass distributions of fragments following de-excitation of 
the giant dipole resonance in 32 S. An illustration of the method's use in collisions is given. 

PACS numbers: 21.10.-k, 21.30.Fe, 21.60.Jz, 27.90.+b 



I. INTRODUCTION 

The nucleus is a complex many-body system which 
shows a rich variety of behaviour, including properties 
that can be identified as single-particle nature as well as 
collective motion. Attempting to describe the nucleus in 
terms of its constituent proton and neutrons and their 
interaction is a difficult task, and one usually focuses on 
a restricted set of degrees of freedom that are of interest 
in the task at hand. 

A venerable approximation that has a fairly general 
application is the self-consistent (or Hartree-Fock) mean- 
field approach, in which the particles are assumed to be 
moving freely (Pauli principle notwithstanding) in a po- 
tential generated by the average effect of the interaction 
of a given nucleon with the rest of the nucleons. A vast 
body of literature exists on approaches starting from the 
mean-field point of view, using a variety of effective nu- 
clear interactions. For recent reviews see e.g. |T] for de- 
tails of mean-field applications of the Skyrme interaction, 
[5] for relativistic models, and [5] for a combined review, 
including the Gogny force. These reviews make clear the 
general applicability of these approaches, and how the 
mean field is a starting point for more elaborate theo- 
ries that bring in correlations beyond the non-interacting 
single-particle picture. 

Standard techniques for going beyond the mean-field 
approach make use of the Hartree-Fock solution as a ref- 
erence state, from which combinations of configuration 
mixing and projection in various ways give rise to corre- 
lated states, and to spectra. Such methods include the 
Random Phase Approximation |U |S] , the generator co- 
ordinate method [51 [7] and shell model-like configuration 
mixing approaches [3 |3] ■ These methods generally rely 
on the outcome of static Hartree-Fock calculations to pro- 
duce the variety of physical states. 

To go beyond the mean-field picture in the time- 
dependent extension of the Hartree-Fock method, a dif- 
ferent approach is usually used. Hartree-Fock can be 
seen, as well as a lowest-order diagram calculation from 
a perturbation theory point of view, as a variational 



principle, and one may instead think about using an 
alternative variational principle to go beyond the limit 
of non-interacting particles. The same is true of Time- 
Dependent Hartree-Fock, whose equations can be derived 
from a variation of the action. In a series of papers, 
Balian and Veneroni developed a general variational prin- 
ciple which included as a special case the time-dependent 
Hartree-Fock approach, but which could be extended in 
various ways [THJ HU HH H3] ■ This form of beyond- mean- 
field approach leads to a natural way to work with time- 
dependent problems. 

Applications are made, in the present work, to cal- 
culations of number fluctuation following decay of giant 
resonance and reaction residues. The method is reviewed 
in secti on |H] some issues of implementation discussed in 
section |III| and results in section |TV] For further details, 
readers are referred to previous work O \TE\ . 



II. THE BALIAN- VENERONI METHOD 

The Time-dependent Hartree-Fock (TDHF) equations 
can be derived (amongst myriad other ways) from the 
variation of the time-dependent action 



5 / (^\id t - H\V) dt = (1) 
Jt ± 

whence the supposition of the state vector being of the 
form of a Slater Determinant containing a set of time- 
dependent single-particle states, {</>}, yields the usual 
TDHF equations 

- £> fc (r,f) f dr / ^(r , > *)t;(r-r , )0 i (r / > t). 
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Choosing a Slater Determinant for the wavefunction, 
or equivalently a one-body density matrix satisfying p 2 = 
p results in an optimal treatment of the dynamics when 
one is concerned with one-body observablcs, with the 
choice of variational principal giving rise to a collective 
path which follows the minimum energy (within the con- 
straints of the wavefunction Ansatz) . There is no reason 
to expect that any other two-body observables be well re- 
produced under these circumstances. If one considers as 
a simple two-body observable the fluctuation in particle 
number 



R c [fm] (N) AiVTDHF A(jV T DHF) m ax AjV B y Change 
8.0 26.637 2.022 2.113 2.351 +16% 

8.5 26.734 2.009 2.098 2.330 +16% 

9.0 26.898 1.988 2.071 2.292 +15% 



TABLE I: Dependence of the bounding radius, R c , separating 
the nucleus from the environment, on the number of nucleons 
in the residual nucleus, (N), the TDHF value of the parti- 
cle number fluctuations, ATVtdhf, the theoretical maximum 
value that TDHF can give Q, A(A^TDHF)max, the Balian- 
Veneroni fluctuation, AA^bv and the change from TDHF to 
BV. 
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one can demonstrate [16] that TDHF can only achieve 
a theoretical upper limit of 
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which is a limit of the method, though the physical mass 
fluctuation may (and does) exceed this. In their work 
[TU] Balian and Veneroni tackled this issue, formulating 
instead a variational principle which sought to include 
the observable of interest in the variational space, so that 
one might trust calculations for at least that observable. 
Their arguments lead to finding stationary values of a 
characteristic function 



J = Tr A(ti)D(ti) 
dt | Tr 



A(t) 



dD(t) 
dt 



-h(A(t),D(t)))(5) 



in which A — exp(— ciV), where N is the one-body 
observable whose square is of interest, whose value at 
t\ is required, D is a density matrix, whose value at 
to is known, e is a small parameter and h is a pseudo- 
Hamiltonian, defined as 



h{A(t),D(t)) = -iTr(A(t) H(t),D(t) 



(6) 



In the special case that one restricts D to be a one- 
body density matrix, one obtains finally an expression 
for the fluctuation of that operator 



is subject to the boundary conditions (|8| at t\. In prac- 
tice, one must perform a TDHF calculation forwards in 
time for the reaction of interest until time i 1; apply the 
transformation (pj) , then evolve a backwards in time to 
t to evaluate (|7^ 



III. IMPLEMENTATION 

The implementation of the Balian-Veneroni (BV) tech- 
nique took as its starting point a three-dimensional time- 
dependent Hartree-Fock code which uses the full Skyrme 
interaction, including time-odd terms and the extended 
spin-orbit force. It has recently been applied to giant 
resonances [TTJ HB] and collisions [T5] [2D]. Since the 
first semi-realistic nuclear TDHF calculations [21], ad- 
vances in computing have recently allowed a more ex- 
tensive and realistic series of calculations [HO [55] [53J [21] . 
Since the BV technique is more computationally intensive 
than TDHF realistic calculations have previously been 
difficult, though previous calculations have been made of 
monopole resonances [25] and reactions [2B], though al- 
ways with some restriction in either spatial symmetry or 
the effective interaction. 

The Balian-Veneroni calculations drive the TDHF 
code by running from to to t\, applying the transfor- 
mation ([8| to the wavefunctions, then running back from 
ti to to. This process is repeated for several values of 
the parameter e in order to take the limit (|7|. The initial 
conditions at to are set up for the collective mode of inter- 
est. In the case of a giant dipole resonance, the starting 
wavefunctions are obtained from the results of a static 
Hartree-Fock calculation, modified by an instantaneous 
boost of the form 



AN 2 (h) = hm ^Tr [(p(t ) - a(t )) 2 } (7) 

where 

<r(ti) = e^p(h)e-^. (8) 

In these equations, p and a are one-body densities 
which evolves according to the TDHF equations a 



<P{t Q ) - e iF *W<f>(to) (9) 

where F T = —1/N for neutron orbitals and 1/Z for 
proton orbitals. f(r) is a spatial formfactor to limit the 
boost to the region of the nucleus. 

Splitting up space into different regions - the nucleus 
and the environment - is key to the Balian-Veneroni pro- 
cess. The wavefunctions are represented in a discretised 
box in space. The box is divided into regions in a way 
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FIG. 1: Logarithmic plot of the total particle density in a head-on collision of two Ca nuclei at a lab energy of 278 MeV. 
Note that the starting position for the collision involves a slight overlap of the nuclei. This is exaggerated by the logarithmic 
scale and the nucleon density at (0, 0) is of order 10~ 5 p c where p c is the density at the centre of one of the nuclei. The decay 
by emission of particles can be seen following the collision. 



specified by the problem at hand. For a single nucleus 
in a giant resonance state, decaying by particle emission, 
a spherical region centered around the nuclear centre of 
mass is considered to be the nucleus, while the rest of the 
box is the environment. This allows for both the num- 
ber of emitted particles to be followed and defined, and 
for the transformation ^ to be applied only to the final 
nucleus and not the emitted flux. Dependence of the re- 
sults of the procedure on the size of the region holding 



the nucleus is one of the several checks for stability that 
has been made. 

Table U shows the results of the variation of the cut- 
off radius for the case of a giant dipolc resonance in 
32 S. Clearly there is some dependence on the cutoff, as 
one would expect. Increasing the size of this radius in- 
creases the number of final particles counted in the nu- 
cleus. However, the changes are small compared with 
the number of dripped particles (i.e. around five) and 
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FIG. 2: Time-dependent dipole moments of S following a 
dipole boost. 

the size of the fluctuations changes rather little. The 
increase in the fluctuation in going from TDHF to BV 
is also extremely consistent. A large amount of testing 
work has checked consistency as a function of geometry 
of cutoff, final running time t\, method of extrapolation 
of the limit e — > and discretisation of the box. For 
further details, see [27] . 

The implementation of collisions relics on bcgining 
with two static HF calculations to give the two collid- 
ing nuclei. These arc placed in the box, and a Galilean 
boost applied to both to set them on a collision path with 
the desired geometry. Following collision, a single frag- 
ment is identified, the BV transformation applied, and 
time run backwards as usual. 

As a final example of an issue of implementation, Fig- 
ure [T] shows a logarithmic plot of the nuclear density 
during the head-on collision of two 40 Ca nuclei. The 
plot shows that following the collision, the excitation of 
the compound nucleus results in the emission of parti- 
cles, seen as the outgoing waves in the third and fourth 
frames. One also sees that at later times the nuclei are 
surrounded by a gas of nucleons formed by reflection of 
the outgoing flux from walls of the box. This unphysi- 
cal reflection is in principle a source of error in the BV 
technique since the interference of the reflected flux with 
the nuclear region can give perturbations to the observ- 
ables. Techniques are available which can eliminate the 
flux [25J|Mj, but they are either computationally costly 
or are not irreversible, or both. 



IV. RESULTS 

Following the prescription of the previous chapter, 
some results of giant dipole resonances in 32 S and col- 
lisions between two 16 nuclei are presented. Figure [2] 



shows the dipole repsonse in 32 S, as calculated with the 
SLy6 Skyrme parameterisation |30j following an isovector 
dipole boost. The moments in the y- and z- directions 
are superimposed since the nucleus is axially symmetric 
with a prolate deformation. The x-axis is the major axis 
and hence the giant resonance has a slower mode in the x 
direction, and two faster modes in the y and z directions, 
as expected from a prolate nucleus. 

Following the deexcitation, the BV transformation is 
applied within the cutoff radius, and the modified density 
matrix evolved back to t^. Results are presented in Table 
[IT] for the fluctuations in particle number for this calcula- 
tion, as well as the equivalent calculation with different 
Skyrme parameterisations. 



Skyrme force 


(N) 


AJVtdhf 


AiV B v 


Change 


SkM* 


26.172 


2.101 


2.287 


+9% 


SLy4 


26.293 


2.060 


2.537 


+23% 


SLy4d 


26.421 


2.053 


2.463 


+20% 


SLy6 


26.637 


2.022 


2.351 


+16% 



TABLE II: Number fluctuation following dipole excition in 
32 S using Skyrme force parameterisations SkM* [3T], SLy4 
[30], SLy4d [32] and SLy6 [30]. 

The variation in predictions across different Skyrme 
forces exceeds that due to any of the numerical uncer- 
tainties discussed briefly in the previous section. This 
is reassuring since it means that this method is sensi- 
tive to properties of the effective interaction, and the 
predictions of fluctuations may be able to be used in 
future to further constrain and select between different 
Skyrme forces, adding to the increasing range of prop- 
erties which are used to judge Skyrme parameterisations 
|33j . It should be pointed out that Skyrme parameterisa- 
tions are fitted essentially to the properties of the ground 
states of doubly-magic nuclei and calculations of, for ex- 
ample, giant resonances in open-shell nuclei are essen- 
tially parameter-free. 

As shown in Figure [l] calculations of collisions have 
also been performed. In this present paper, no analysis 
is presented beyond the discussion of the figure in the 
previous section. Work on completing the analysis is in 
progress, and will be reported elsewhere. The first anal- 
ysis of collisions between 16 nuclei can be found in the 
literature [T5] . 

V. CONCLUSION 

The Balian-Veneroni method for reliably evaluating 
the fluctuations of one-body observables has been im- 
plemented for symmetry-unrestricted nuclear dynamics 
with the Skyrme interaction. First results have been pre- 
sented for giant resonances and collisions. The approach 
is computationally intensive and has some inherent nu- 
merical issues involved in its implementation, but these 
are well in control, and the differences evident in different 
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Skyrme parameterisations are larger than the numerical 
uncertainty. The Balian-Veneroni technique can there- 
fore be used both to calculate fluctuations, as well as to 
serve as a filter for Skyrme interactions. 
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